Two-bath model for activated surface diffusion of interacting adsorbates 
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The difTusion and low vibrational motions of adsorbates on surfaces can be well described by a 
purely stochastic model, the so-called interacting single adsorbate model, for low-moderate coverages 
{9 < 0.12). Within this model, the effects of thermal surface phonons and adsorbate-adsorbate 
collisions are accounted for by two uncorrelated noise functions which arise in a natural way from 
a two-bath model based on a generalization of the one-bath Caldeira-Leggett Hamiltonian. As an 
illustration, the model is applied to the diffusion of Na atoms on a Cu(OOl) surface with different 
coverages. 

PACS numbers: 05.10.Gg, 05.40.-a, 68.43.-h 
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I. INTRODUCTION 

Nowadays, quasielastic Helium atom scatteringiiS, 
(QHAS), '^He spin-echo2. and neutron spin-echof^ are 
well established experimental techniques to study fast 
diffusion processes which are aimed at determining gas- 
surface interaction potentials. The observables that can 
be measured with these techniques are the so-called dy- 
namic structure factor and the intermediate scattering 
function (or polarization function), which is the inverse 
time Fourier transform of the former. The diifuse elas- 
tic intensity of the He atoms scattered away at large an- 
gles from the specular direction provides very detailed in- 
formation about the mobility of adsorbates on surfaces. 
Based on the transition matrix formalism, and applied 
to the first technique, Manson and Celli^ proposed a 
quantum diffuse inelastic theory for small coverages of 
adsorbates on the surface, ignoring multiple scattering 
effects with the incoming He atoms. From this theory, 
they obtained the dynamical structure factor under the 
assumption that all crystal vibrational modes and point- 
like scattering centers satisfy the harmonic approxima- 
tion with a given frequency distribution. 

Alternatively, the study and analysis of surface dif- 
fusion processes can also be tackled by means of the 
Langevin formalism in two dimensions. This assumption 
is reasonable because, usually, the adparticle motion nor- 
mal to the surface involves vibrational modes with much 
higher frequencies than in-plane vibrations. Nonetheless, 
some recent work showa^ that motion perpendicular to 
the surface, related to a translational hopping diffusion 
process, could be important for coverages 9 > 0.05, as 
noticed for small parallel momentum transfers. As is 
well known, one of the most popular and ubiquitous phe- 
nomcnological equations for systems interacting with en- 
vironments is the generalized Langevin equation (GLE). 
In the case of surface diffusion, this equation reads as 



where R = (x, y) is the adsorbate position at a certain 
equilibrium distance upon the surface and F = —WV 
is the deterministic force acting on the adsorbate, with 
V being the periodic adsorbate-surface interaction po- 
tential determined from He scattering in the zero cover- 
age and zero surface temperature limits.^ Equation ([1]) 
is a stochastic integro-differential equation with addi- 
tive fluctuations and linear dissipation, for neither the 
fluctuating force, 6Ff = Ff — (F/), nor the dissipative 
kernel, K(i — t'), depend on the adsorbate position R. 
Usually, the fluctuating force is considered as a zero- 
centered Gaussian, completely specified by its autocor- 
relation function, C{t — r), which is assumed to be sta- 
tionary. Because of the fluctuation-dissipation theorem, 
which assures the system reaches an equilibrium asymp- 
totically in time (steady state), we have 



C(i-r) = kBTK{t-T), 



(2) 



R(i) 



K{t - t') R(t') dt' + F[R(t)] + 5F}{t), (1) 



where /cs is the Boltzmann constant and T is the surface 
temperature. Whenever this description is valid, Eq. ^ 
is independent of the particular details of the interac- 
tion between the system and the surrounding heat bath. 
For Ohmic friction (no memory effects), the kernel can 
be expressed as a delta function and Eq. ([T]) becomes a 
standard Langevin equation. 

Former experimental and theoretical studies of Na dif- 
fusion on a Cu(OOl) surface at low coverages were car- 
ried out by assuming that the corresponding dynam- 
ics is well described by the motion of a single adsor- 
bate, governed by the standard Langevin equation. 
The thermal effects induced by the surface temperature 
were accounted for by a Gaussian white noise, this being 
in agreement with contemporary full molecular dynam- 
ics simulations,™ where the motion of both the surface 
atoms and a sample of adatoms scattered over the sur- 
face was taken into account. As coverage increases, the 
Langevin treatment has still been used in the literature;^ 
but in combination with molecular dynamics (Langevin 
Molecular Dynamics, LMD), where the inter-adsorbate 
forces are commonly described by repulsive dipole-dipole 
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potentialsi"°''""' An interesting result obtained from this 
type of simulations is the appearance of ordered struc- 
tures at 6* = 0.125 and — 0.16, compatible with a dom- 
inant repulsive interaction, as observed experimentally..^^ 
Therefore, this fact should be taken into account as a lim- 
itation of the range of validity for fully based stochastic 
approaches, since they can not describe the appearance of 
this kind of collective phenomena or phases. Within the 
context of the GLE, the memory effects on the frictional 
damping in surface diffusion have also been discussecU^ 
by considering a Hamiltonian similar to the Caldeira- 
Leggett Hamiltonianfi^"— This leads to an analysis in 
terms of the relationship between the surface excitation 
frequency (i.e., the Debye frequency), lod, and the char- 
acteristic vibrational frequency of the adatoms, wq. In 
general, when ud is greater than wq, memory effects 
are not very important and, therefore, an Ohmic fric- 
tion can be assumed. This is the case, for example, for 
Na atoms diffusing on a Cu(OOl) surface, where the De- 
bye frequency is about 30 meV, while the lowest T-mode 
frequency is about 6 meV. 

Though the GLE is a phenomenological equa- 
tion, it can be formally derived from a microscopic 
Hamiltoniar^iirJ^ consisting of the system coupled to a 
bath of harmonic oscillators (the one-bath model). This 
formalism has been used, for example, to study quan- 
tum tunneling in dissipative systems j^iiiii^ the so-called 
Kramers' turnover problemji^"^^ and vibrational dephas- 
ing rates 1^ A remarkable aspects of this formalism is 
that it allows us to establish straightforwardly a connec- 
tion between the Langevin formalism in surface diffusion 
and quantum mechanical approaches, such as the one in- 
troduced by Manson and Celh^^ Hence, since Eq. ([1]) is 
very general, we can also assume that there are two in- 
dependent, uncorrelated baths (one for the phonons and 
another for the adsorbates) in order to describe surface 
diffusion processes with interacting adsorbates within the 
framework of a two-bath model based on a generalization 
of the one-bath Caldeira-Legget Hamiltonian. In this 
way, two independent fluctuating forces are present and 
the total kernel then consists of the sum of the kernels 
associated with each fluctuating force. As mentioned, if 
each bath is assumed as Ohmic, both kernels will then 
be describable in terms of delta functions in time. This 
leads immediately to the so-called interacting single ad- 
sorbate (ISA) model^"— which has been shown to re- 
produce fairly well experimental surface diffusion data 
under the presence of interacting adsorbates up to mod- 
erate coverages^ [0 ~ 0.12). Kramers' theory can also 
be generalized using two baths to infer physical prop- 
erties of the diffusing particle through the quasielastic 
(Q) peak obtained from experimental data, as previ- 
ously proposed within the one-bath modeL ^^'^^ Further- 
more, the vibrational dephasing theory proposed for the 
low-frequency frustrated translational motion or T-mode 
within the one-bath model can also be straightforwardly 
extended^^S This allows us to provide some analytical 
expressions, flrmly based on formal grounds, for a fitting 



procedure in order to extract relevant information from 
the experimental data. 

The organization of this work is as follows. In Sec. [H] 
we present the basic building blocks of the theoretical for- 
malism which allow to connect the two-bath and the ISA 
models, as well as their connection to the experiment. In 
order to illustrate the applicability of this formalism, in 
Sec. lIIII we present an analysis of Na diffusion on Cu(OOl) 
where a lot of experimental data are available from the 
literature Finally, in Sec. lIVI we summarize the main 
conclusions derived from this work. 



II. THEORY 
A. Basic formalism 

In diffusion experiments carried out by means of 
QHAS, one measures the differential reflection coefficient 
which, in analogy to liquids^^ reads as 



^^27^(AK,w) 



= ridF jj G(R,t)e*('^^-^-"*)dRdt. 



(3) 



This expression gives the probability that the probe He 
atoms scattered from the diffusing collective (distributed 
upon the surface with a concentration n^) reach a certain 
solid angle with an energy exchange hjj = Ef — Ei 
and wave vector transfer parallel to the surface AK = 
— K^. In Eq. F is the atomic form factor, which 
depends on the interaction potential between the probe 
atoms in the beam and the adparticles on the surface. 
On the other hand, S'(AK,a;) is the so-called dynamic 
structure factor or scattering law. Within the context of 
the linear response theory, this observable corresponds 
to a response function and provides a complete infor- 
mation about the dynamics and structure of the adsor- 
bates through particle distribution functions. Experi- 
mental information about long-distance and long-time 
correlations are obtained for small values of AK and 
small energy transfers huj, respectively, through the Q 
and T-mode peaks. Particle distribution functions can 
be well-described by means of the so-called van Hove 
or time-dependent pair correlation function,-?^ G{Il,t), 
when dealing with interacting particles. Given a particle 
at the origin at some arbitrary initial time t = 0, G(R, t) 
gives the averaged probability of finding the same or an- 
other particle at the surface position R at a time t. This 
function thus generalizes the well-known pair distribution 
function ^(R), commonly used within the context of sta- 
tistical mechanics by providing information about 
the interacting particle dynamics. 

Alternatively, the dynamic structure factor can be 
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expressedaii as 

= J e-"^' I{AK,t)dt, (4) 

where the brackets in the integral denote an ensemble 
average. In the second line of Eq. Q, 

I{AK,t) = (e"*'^K-[^(*)-^(o)l) = ^g-iAK-/o%AK(t') dt'^ 

(5) 

is the intermediate scattering function (also known as 
polarization^^— in spin-echo experiments), which is the 
space Fourier transform of G{'R.,t). In Eq. (I5|), wak is 
the velocity of the adparticle projected onto the direction 
of the parallel momentum transfer, AK. The averages 
involved in Eq. ([5]) can be obtained from Langevin nu- 
merical simulations for an adparticle in the presence of 
an external field of force V^(R), which gives the substrate- 
adsorbate interaction, and two uncorrelated noise func- 
tions leading to two friction coefhcients: one associated 
with thermal noise (surface phonons) and another with 
the collisional friction due to the collisions among adsor- 
bates. 



B. Kramers' theory for interacting adsorbates 

As mentioned above, Manson and Celli^ proposed a 
quantum diffuse inelastic theory for small and interme- 
diate coverages based on the transition matrix formal- 
ism, ignoring multiple scattering effects with He atoms. 
Within this approach, 5(AK, cj) is obtained after assum- 
ing that all crystal vibrational modes {N) and point-like 
scattering centers (M) satisfy the harmonic approxima- 
tion with a given frequency distribution. In a similar way, 
we can describe the diffusion of interacting adsorbates by 
means of two types of independent, uncorrelated baths. 

Based on Kramers' theory ; '^^i^^ a quantum and classi- 
cal theory of surface diffusion at very low coverages has 
been developed in recent years i^^i^iSl Within this the- 
ory, the corresponding GLE is split up into two coupled 
equations, one for each degree of freedom. This GLE 
arises from a total system-|-bath Hamiltonian expressed 
in terms of a single bath consisting of a set of N harmonic 
oscillators. This bath simulates the surface thermal ef- 
fects on the adsorbate. For low and moderate coverage 
{9 < 0.12), a similar model Hamiltonian but with two 
baths can also be formulated. In this case, apart from 
the surface thermal effects, we take one adsorbate as the 
tagged particle or system, while the remaining ones con- 
stitute the second bath, which is described as a collec- 
tion of M harmonic oscillators. Obviously, the charac- 
teristics of the second bath are going to be dependent on 
the surface coverage considered in the experiment. Thus, 
the corresponding total (system-f two-bath) Hamiltonian 



reads as 
H = 



pi+Pv 



2m 
1 ^ 

1 ^ 

4e 



niiLUxi 



i=l 
M 



1 



p^- 



niiUJy, 

dr.. 



m-j (jJx-Xj 



mjUJxj 



, (6) 



where V{x, y) is the periodic potential describing the sub- 
strate and [pxiPy) and {x,y) are the adparticle momenta 
and positions on the surface, and {px^^Xi) and {pyi,yi) 
with i = I,-- - ,iV are the momenta and positions of 
the ith bath oscillator (phonon), with masses and fre- 
quencies given by mi and uik^ {k = x,y), respectively. 
Phonons with polarization in the z-direction are not con- 
sidered. The same holds for the barred magnitudes, 
though they are associated with an assumed harmonic 
oscillator bath of M adsorbates, thus neglecting eventual 
long-distance correlations. The Cfe^ and dk^ coefficients in 
both directions {k = x, y) give the adsorbate-substrate 
and adsorbate-adsorbate coupling strengths, respectively. 

The spectral density for the two baths is defined simi- 
larly to the case of one bath. 



N 



O 2 ^(^-WfcJ 

2 ^ m.ujl 



M 

-^V - 

2 rujujl 



6{uj-ujkj), k^x,y,(7) 



though now it contains two terms: one spectral density 
is associated with the surface phonons and the other one 
with the bath of adsorbates. Equation ([7]) enables pass- 
ing to a continuum model provided the time-dependent 
friction can be expressed as 



mit) 



4- ■^ki^'l ^ 1 IQ\ 

coswr dw, k — x,y. (8) 



In the case of Ohmic friction, i.e., ri{t) = 2ri6{t), Eq. ([T]) 
reduces to the standard Langevin equation^"— (the 
delta function counts only one half when the integration 
is carried out from to oo). 



R 



-?7R + F(R) 



(9) 



which constitutes the basis of the ISA model. In this 
equation, SF f is given by the sum of two noncorre- 
lated noises: the lattice (thermal) vibrational effects 
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and the adsorbate-adsorbate collisions, which are sim- 
ulated by a Gaussian white noise (G) and a white shot 
noise (S), respectively. Thus, for each degree of free- 
dom, we have SFf{t) = 5Ff {t) + 5Ff{t). The Gaus- 
sian white noise satisfies the properties: {F'^(t)) = and 
{Ff{t)Ff{t')) = 2mjkBTS{t' -t), where m is the adsor- 
bate mass and 7 is the (constant) friction coefficient mea- 
suring the adsorbate-phonon coupling strength. On the 
other hand, the shot noise is given'^* by a sum of pulses 
mimicking the collision impacts. In the Markovian ap- 
proximation, these collisions are assumed to be sudden 
(strong but elastic), with post-collision effects relaxing 
exponentially at a constant rate much larger than the 
average number of collision per time unit or collisional 
friction, A. The memory function or kernel associated 
with the shot noise in ([T|) then becomes local in time and 
6Fj{t) will display features of a white shot noise. In 
other words, the adsorbate-adsorbate collisions are de- 
scribed by means of a white shot noise as a limiting case 
of a colored shot noise. Thus, K{t — t') ~ 2rj S{t — t'), 
where the total (Ohmic) friction is 77 = 7 + A. A simple 
relationship between the collisional friction. A, and the 
coverage, 9, at a temperature T is givei*^ by 



A = 



(10) 



The theory of activated surface diffusion in one di- 
mension was developed from Kramers' solution32i2I to 
the problem of escape from a metastable well,^i2& The 
underlying dynamics is assumed well described by the 
Langevin equation provided that the reduced barrier 
height is of the order of ~3 or higher (i.e., /UbT » 3). 
The energy loss to the bath of trajectories close to the 
barrier top is given by classical mechanics and the po- 
tential at the barrier top is approximately parabolic. 
Kramers' based theory with finite barrier correction 
terms can then be replaced by Langevin numerical sim- 
ulations. As is well known, the so-called turnover 
regiorM^ is observed when the transmission factor or 
prefactor of the exponential law for the escape rate is 
plotted versus the friction coefficient at a given surface 
temperature. Such a factor behaves linearly with 77 for 
relatively low frictions (i.e., ~ 77) and goes like ~ 77"^ 
at the low-to-moderate friction regime (Smoluchowski 
limit). In between, it passes through a maximum, thus 
undergoing the turnover, where it displays a characteris- 
tic smooth shape. 

The starting point of the Kramers' model is a kinetic 
equation in one dimension for the stationary flux of par- 
ticles exiting each well at either barrier. This flux is 
affected by the rate of particles exiting the jth well and 
those arriving at the well from the two neighboring wells, 
j + 1 and j — 1 . Here we are going to give the main analyt- 
ical expressions derived from Kramers' theory, but more 
details can be found elsewhere i^^i^^'^^ A central quantity 
in the theory is the reduced average energy loss, 5, to 
the bath as the adatom traverses from one barrier to the 



next. If a single cosine potential with a barrier height 
— 2Vo is considered, the energy loss is given by 
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8V0V 
ksTujo ' 



(11) 



where ojq — 2Try^VQ /ma^ is the harmonic frequency of 
the oscillations near the well bottom, m is the mass of 
the adatom and a is the unit cell length. Since many 
experiments or calculations are carried out typically un- 
der conditions of large reduced barrier heights, S can be 
unity or even larger, even though the damping constant 
is rather small. 

In the moderate to strong friction regime, where the 
rate limiting step is the spatial diffusion (sd) across the 
barrier, the rate of the escape from the well in both direc- 
tions is given by the Kramers-Grote-Hynes formula^^i^S 



A* 



LJ* IT 

where the corresponding prefactor is 



1 



(12) 



(13) 



the magnitudes associated with the barrier being denoted 
by f . This prefactor appears as a renormalization taking 
into account recrossings, since we are working implicitly 
in normal mode coordinates for the diffusing particle and 
the baths. Finally, for the partial rates one finds 




cos{jAK) 



1 + P^(x) - 2P{x) cos(Aii') 



dx 



(14) 



The escape rate from the ze- 



X dAK, 

where P{x) = e-*/4™s^(^). 
roth surface well is k = — Fg and the relative probability 
for a jump of length j is given by the probability of trap- 
ping at the jth well, Pj — Tj/n. For a one-dimensional 
periodic potential, the diffusion coefficient is related to 
the escape rate as 



2 ^ ' 2 



(15) 



where (P) is the mean square path length. By using 
Eq. this diffusion coefficient can be reexpressed in 

a more compact form as 



exp 



^ ' ln[l + P(a;)]da;| 



(16) 



where D^d — (l/2)a^Fsc( is the diffusion coefficient in the 
spatial diffusion regime and T is the depopulation factor 
for the metastable well. 



T — exp 



l^^'^'ln [1 - P(x)] dxj , 



(17) 
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first given by Melnikov.-2fi 

Now, in analogy to other models, such as the Chudley- 
EUiott niodelf2S,d2 an analytical expression for the full 
width at half maximum (FWHM) of the dynamic struc- 
ture factor can also be obtained by imposing a master 
equation for the intermediate scattering function<^ If 
the dynamic structure factor has a Lorentzian shape, the 
FWHM is given by 



T{AK) = ATsd sin' 



X exp < — 



7r/2 



In 



l + P'^{x)-2P{x) cos(Aif) 



dx 



(18) 



This equation is important in the sense that, assuming 
the validity of Kramers' model and the master equation 
approach mentioned above, it allows for a direct com- 
parison with the experimental data and/or Langevin nu- 
merical simulations and therefore an estimation of the 
spatial diffusion rate and the energy loss d. From 
these parameters and their temperature dependence, one 
can further infer the barrier height, the friction coefhcient 
and the barrier frequency. 




FIG. 1: FWHM as a function of the parallel momentum trans- 
fer for Na diffusion on Cu(OOl) along the [100] direction at 
T = 300 K. Experimental dataii are indicated with symbols, 
while the solid lines are the fitted curves, obtained by applying 
Eq. (|18|l . Two different values of the coverage are considered: 
01 = 0.028 (black) and 62 = 0.078 (red). 



location, (ujt), is given by 

{ujt} = t^cff + ksT 1 - ksT — 3- 



(20) 



C. Vibrational dephasing and T-mode 

The vibrational dephasing of an adsorbate on a surface 
arises from the random frequency fluctuation undergone 
by the presence of the surface bath. Dephasing mainly 
emerges from the anharmonicity of the potential describ- 
ing the adsorbate vibrations. When the influence of the 
coverage becomes important, another source of fluctu- 
ations arises, which is associated with the interaction 
among adsorbates. Within the ISA model^^"— for ex- 
ample, this effect is accounted for by a white shot noise. 
Dephasing times or rates can be obtained from the auto- 
correlation function of the frequency fluctuation through 
the Green-Kubo expression.™ 

In Ref. HO, simple analytic expressions are derived for 
the vibrational line shapes of the T-mode from a Hamil- 
tonian analogous to ([6]), though in the case of the one- 
bath model. In particular, these expressions describe the 
location and FWHM of the T-mode peak as a function 
of the surface temperature. According to such a model, 
to leading order, the T-mode potential of mean force can 
be expressed generically as 



V{q) = i Lulq^ + IUq\ 



(19) 



where q is the mass weighted T-mode vibrational coordi- 
nate, ujQ is the harmonic frequency, and if 4 is the anha- 
monicity constant. To second order in the anharmonic- 
ity and provided the damping constant 77 is not too large 
iv/^o ^ 2), for Ohmic friction one finds that the peak 



where lucS = \/oJq — {r//2)^ in the case of the two-bath 
model described by ([6]). Similarly, the T-mode peak 
FWHM, CT, depends on temperature and friction as 



G\Ki\ 

OJqUJcS ' 



(21) 



III. RESULTS 

Kramers' turnover theory provides a two-parameter 
description of the diffusion process through the FWHM 
of the Q-peak versus the parallel momentum transfer. 
As seen in the preceding Section, these two parameters 
are the spatial diffusion rate and the energy loss. This 
treatment suits particularly well in fitting procedures ap- 
plied to experimental data and/or numerical Langevin- 
type results due to the closed analytical expressions that 
are involved. Very good estimates of rates, diffusion co- 
efficients and jump distributions are obtained when the 
barrier for diffusion exceeds the thermal energy and fi- 
nite barrier corrections are also included. In Fig. [l] 
experimentalii (symbols) and numerically fitted (solid 
lines) FWHM are plotted as a function of the parallel mo- 
mentum transfer AK for Na diffusion on Cu(OOl) along 
the [100] direction at T = 300 K. Two different cov- 
erages are considered: 9i = 0.028 (black squares/line) 
and 6*2 — 0.078 (red circles/line). Fittings have been 
carried out using Eq. (IT5|) . The fitted values obtained 
for the spacial diffusion rate are r^^j = 0.250 meV and 
0.238 meV, and for the energy loss, 5i = 1.29 and 
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FIG. 2: T-mode frequency as a function of the surface tem- 
perature for Na diffusion on Cu(OOl) along the direction 
[100]. Experimental datai^ are indicated with symbols, while 
the solid lines are the fitted curves, obtained by applying 
Eq. (|20p . Three different values of the coverage are consid- 
ered: Oi = 0.028 (black), 02 = 0.078 (red) and 63 = 0.125 
(blue). 



$2 = 1.97. The broadening of the Q-peak with cover- 
age is quite well reproduced by Eq. Even more, as 
is clearly seen from Eq. (IT6l) . the diffusion coefficient de- 
creases with the coverage (as shown previously in Ref.[26l) 
and the energy loss increases with the coverage or friction 
according to Eq. ([TT|) . 

The T-mode peak position and its FWHM are plotted 
in Figs. [2] and [SI respectively, as a function of the surface 
temperature and for three coverages: 9i = 0.028 (black), 
6*2 = 0.078 (red) and 63 = 0.125 (blue). In Fig. [2 the 
best fitted curves (solid lines) with respect to the corre- 
sponding experimental data"^^ (symbols) for the T-mode 
peak positions have been obtained from Eq. (|20p . The 
corresponding frequencies obtained when Eq. (1201) is ap- 
plied are: w^^^ = 2.13 x 10"", w^^^ = 2.3 x 10"*, and 

ojg'^'' = 2.4 X 10^'', all given in atomic units (the nom- 
inal value being 2.20 x 10^^ a.u.). With the coverage, 
the T-mode frequency shifts towards higher values. Af- 
ter the definition of the effective frequency in Eq. ([20l) , it 
should display the opposite trend with the coverage since 
a higher friction coefficient is required in the ISA model. 
However, the changes on ujq are so small (less than 10%) 
that we can not rule out small changes in the mean force 
potential. On the other hand, the best fitted values found 
for the anharmonicity constants (also through Eq. (|20)) ) 

are ^ = -6.17 x lO^^^, ^ = -8.05 x lO^i"*, and 

K^f ^ = —9.25 X 10"^"*, also in atomic units. For a co- 
sine potential, is negative and therefore the anhar- 
monicity produces a red shift with a linear temperature 
dependence.-^-fi The agreement with the surface tempera- 
ture is fairly good when second-order corrections in the 
anharmonicity are included. The different loq and K4 val- 
ues are slightly modified due to coverage, thus indicating 
a modification of the barrier/ well and anharmonicity due 
to increasing amount of adsorbates, respectively. In gen- 
eral, in the presence of only one bath, the new system 
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FIG. 3: T-mode FWHM as a function of the surface tempera- 
ture for Na diffusion on Cu(OOl) along the direction [100] . The 
solid lines have been obtained by means of Eq. pif) and using 
the values arisen from the fittings displayed in Fig. [21 Three 
different values of the coverage are considered: 61 = 0.028 
(black), 62 = 0.078 (red) and 63 = 0.125 (blue). 

plus bath frequencies are renormalized due to the sur- 
face friction or coupling with the substratci^ This renor- 
malization takes into account the sum over all harmonic 
oscillators of the bath. In the two bath scenario, the 
corresponding renormalization has an extra sum coming 
from the adsorbates. Finally, Fig. [3| remains as a predic- 
tion since no experimental information is available about 
FWHM of the T-mode peaks. However, what is well- 
known from experiments is that the T-mode peak broad- 
ens with the coverage, which is theoretically predicted by 
Eq. 

IV. CONCLUSIONS 

The main physical consequences of the approach pre- 
sented here are that, on average, the tagged adpar- 
ticle is seen as a diffusing particle among adsorbates 
which modify effective barriers and wells leading to the 
renormalization of frequencies as showed from a nor- 
mal mode analysis.— The corresponding renormalized 
frequencies have been reported for the one-bath model 
elsewhere»i^ A straightforward generalization of the effec- 
tive frequencies can be done for the two-bath model. In 
particular, this new approach could provide a theoretical 
framework to better interpret the coverage dependence 
of tunneling diffusion experiments^ within a quantum 
Markovian formalism.— Many times, in surface diffusion 
LMD simulations surface barriers are only changed with 
the coverage in order to reproduce experimental data, 
keeping the corresponding frequencies fixed. 
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